Exact asymptotics of monomer-dimer model on rectangular semi-infinite lattices 



Yong Kon£| 
Department of Mathematics 
National University of Singapore 
Singapore 1 1 754 3 

(Dated: February 6, 2008) 

By using the asymptotic theory of Pemantle and Wilson, exact asymptotic expansions of the free 
energy of the monomer-dimer model on rectangular tixoo lattices in terms of dimer density are 
obtained for small values of n, at both high and low dimer density limits. In the high dimer density 
limit, the theoretical results confirm the dependence of the free energy on the parity of n, a result 
obtained previously by computational methods. In the low dimer density limit, the free energy on 
a cylinder n x oo lattice strip has exactly the same first n terms in the series expansion as that of 
infinite oo x oo lattice. 
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I. INTRODUCTION 

In the previous investigation of the monomer-dimer 
model on two-dimensional rectangular lattices by com- 
putational methods, it is found that at high dimer den- 
sity, the free energy depends on the parity of the width of 
the lattice strip (for a review of monomer-dimer model, 
see Ref . Q], H and references cited therein) . For a m x n 
lattice, we define the dimer density p as p = 2s /mn, 
where s is the number of dimers. With this definition, a 
close-packed lattice will have p = 1. The grand canoni- 
cal partition function of the monomer-dimer system in a 
m x n two-dimensional lattice is 
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where a s (m, n) is the number of distinct ways to arrange 
s dimers on the m x n lattice, a„ hn (p) is the corresponding 
number at the fixed dimer density p, x is the activity of 
the dimer, and N is the maximum possible dimer on the 
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lattice: N = \mn/2\ . The free energy per lattice site at 
a given dimer density p is defined as 
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fm,n(p) 



For a semi-infinite ooxn lattice at high dimer density 
limit (p — > 1), the computational method shows that the 
free energy foc.n(p) depends on the parity of the lattice 
width n [2j: 



/oc,n(p) ~ Z^" iCC (l) - 



(1 — p) ln(l — p) n is odd 
, |(1 — p) ln(l — p) n is even 

(2) 

where /^" lco (l) is the free energy of close- packed lattice 
with width n, the exact expression of which is known 
and is dependent on the boundary conditions of the lat- 
tice (cylinder or with free edges) [1, 0| • A slightly more 
general expression, for a finite mxn lattice at high dimer 
density limit as p — > 1 and m — * oo, is given by |2| 
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For the coefficient of the term (1 — p) ln(l — p) in the 
above expansion of the free energy foo.n(p), computa- 
tional results give values very close to the value —1 when 
n is odd, and —1/2 when n is even. Although these nu- 
merical evidences strongly support the coefficients of the 
parity-depending term (1 — p) ln(l — p) as shown in Eq. 
(J3j> , there still exists slight possibility that the observa- 
tions be explained in another way. This is mainly due to 
the following two reasons. Firstly, the sequence of n used 



in the computational studies is not very long. For cylin- 
der lattices, the largest n used is 17. For lattices with 
free edges, maximum value of n — 16 is used. Secondly, 
at high dimer density limit, the convergent rate is the 
poorest and some heuristic weighting averages had to be 
used in the fitting procedure L 2j . Due to these uncertain- 
ties, it might be possible that the different coefficients 
observed for these small number of even and odd lattice 
widths are just the initial values of two slowly changing 
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sequences. When the lattice width n is small, the devia- 
tion of the coefficient from —1 or —1/2 is so small that it 
could not be detected numerically. As n becomes bigger, 
these two sequences might converge to some intermediate 
values between —1 and —1/2. 

In this report, we use the asymptotic theory of Pe- 
mantle and Wilson Q to get exact asymptotics for the 
free energy foo.n(p) of n x oo lattices. The exact asymp- 
totics show that the coefficients of (1 — p)ln(l — p) in 
the high dimer density expansion of /oo,n(p) is exactly 
— 1 when n is odd, and —1/2 when n is even. Further- 
more, the Pemantle and Wilson method is also used to 
investigate the low dimer density expansion of foo.n(p)- 
In the low dimer density limit, the free energy on a cylin- 
der n x oo lattice has exactly the same first n terms as 
that of infinite oo x oo lattice. In addition, closed form 
expressions can be obtained for the differences of the co- 
efficients between finite and infinite lattices of the next 
two terms [Eq. (p~5|) and Table IJ. These properties not 
only explain the fast convergence rate of the free energy 
on cylinder lattices when dimer density is low, but also 
provide a quantitative indicator of the errors when the 
results of finite size lattices are used to approximate the 
infinite lattice. 

The article is organized as follows. In Section [III the 
Pemantle and Wilson (PW) method for asymptotics of 
multivariate generating functions is summarized. The 
starting point for the PW method is the multivariate 
generating function of the model under study, and in our 
case of the monomer-dimer model, the generating func- 
tions are bivariate. In Appendix, the bivariate generating 
functions of monomer-dimer models in two-dimensional 
rectangular lattices are listed for small values of n. These 
generating functions are used as the input to the PW 
method in this article. In Section HTT1 the exact asymp- 
totic expansions of foc.n(p) at high dimer density is de- 
rived for some small values of n. The coefficients obtained 
for the (1 — p) ln(l — p) term confirm the dependence on 
the parity of n, as shown in Eq. ([3]). In Section [TVl we 
discuss the exact asymptotic expansions of foc,n{p) at 
low dimer density. 



II. PEMANTLE AND WILSON METHOD FOR 
ASYMPTOTICS OF MULTIVARIATE 
GENERATING FUNCTIONS 

To extract asymptotics from a sequence, it's usually 
useful to utilize its associated generating function. The 
method for extract asymptotics from univariate gener- 
ating function is well-known For multivariate gen- 
erating functions, however, the techniques were "almost 
entirely missing" until the recent development of the Pe- 
mantle and Wilson method For combinatorial prob- 
lems with known generating functions, the method can 
be applied in a systematic way. The theory applies to 
generating functions with multiple variables, and for the 
bivariate case that we are interested here, the generating 



function of two variables takes the form 
F(x,y) 



G(x,y) 



H(x,y) 
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where F(x, y) and H {x, y) are analytic, and H (0, 0) ^ 0. 
In this case, PW method gives the asymptotic expression 
as 
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where (xq, yo) is the positive solution to the two equations 

rr/ . dH dH 

H{x,y)=0, mx — = ( 6 ) 

and Q(x,y) is defined as 

Q(x,y) = - (xH x )(yH y ) 2 - (yH y )(xH x ) 2 

- [(yH y ) 2 {x 2 H xx ) + (xH x ) 2 (y 2 H yy ) 
-2{xH x ){yH y ){xyH xy )}. 

Here H x , H y , etc. are partial derivatives dH/dx, dH/dy, 
and so on. One of the advantages of the method over pre- 
vious ones is that the convergence of Eq. ([5]) is uniform 
when s/m and m/s are bounded. 

For the monomer-dimer model discussed here, with n 
as the finite width of the lattice strip, m as the length, 
and s as the number of dimers, we can construct the 
bivariate generating function G n (x,y) as 



G n (x,y) 
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For the monomer-dimer model, as well as a large class of 
lattice models in statistical physics, the bivariate generat- 
ing function G(x, y) is always in the form of Eq. ([4]), with 
F(x,y) and H(x,y) as polynomials in x and y. In fact, 
for monomer-dimer model as well as other lattice mod- 
els, we can get H(x,y) directly from matrix M n used in 
the recursive formula to calculate the partition functions 

Here the vector Q m consists of the partition function Eq. 
([T]) as well as other contracted partition functions Q. 
The function H(x, y) is closely related to the character- 
istic function of M n Q: H(x, y) = det(M„ - I/y) x y w , 
where w is the size of the matrix M n . More discussions 
can be found in the Appendix. 

The relation between the number of dimers s and the 
dimer density p is given by s = pmn/2. If we fix the 
dimer density p, and substitute this relation into Eq. ((SJ), 
then we see that the solution (xo, yo) of Eq. © depends 
only on p and n, and does not depend on mors: 



H{x,y) = 0, 

dH pny dH 



dx 



2 dy 



(8) 
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Substituting this solution (xo(n, p), yo(n, p)) into Eq. (O 
we obtain 

, / \ 1 w pn/2 , Hum 

fm,n(P) ln ( x VO - o 

n 2 ran 
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From this asymptotic expansion we obtain the logarith- 
mic correction term with coefficient of exactly —1/2 (the 
second term in Eq. ©), for both even and odd val- 
ues of n. In fact, the PW asymptotic theory predicts 
that there exists such a logarithmic correction term with 
a coefficient of —1/2 for a large class of lattice models 
when the two variables involved are proportional, that 
is, when the models are at a fixed "density" . For those 
lattice models which can be described by bivariate gen- 
erating functions, this logarithmic correction term with 
a coefficient of —1/2 is universal when those models are 
at fixed "density" . For the monomer-dimer model, this 
proportional relation is for s and m with s = pmn/2. 

When the dimer density p and the lattice width n are 
fixed, the first term of Eq. §§§ is a constant and does not 
depend on m. We identify it as foo,n(p) 



foo.nip) 
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In theory, as long as the bivariate generating function 
G(x,y) or its denominator H(x,y) is known, (xq, yo) 
could be solved from the system of equations Eq. ([5]) and 



foo,n(p) could be obtained from Eq. (fTOjl . In practice, 
however, only very small values of n can be handled this 
way. When n = 1, xq and yo can be solved as 



x 



which leads to 
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This expression is exact for < p < 1. 

When n = 2, for both the cylinder lattice and the lat- 
tice with free boundaries, H(x,y) is a cubic polynomial 
(Appendix). From Eq. ([8]) a quartic equation can be 
obtained for xq. For the cylinder lattice, xq satisfies the 
following quartic equation 



32 (l-p) 3 a; 4 + 144 (1 - p)° x 

+ 4 (l-p)(l0p 2 -20p + 3) x 2 
+ 4 (2-p) (3p 2 -3p+l)x 

-p(l-p)(2-p) = 0. (11) 

After Xq is solved, yo can be obtained as a rational func- 
tion in p as 
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where 



w(x ,p) = (1440 p 5 - 11424 p 2 + 12960 p 3 - 672 + 4672 p- 6976 p 4 ) x 3 

+ (-3360 + 60432 p 3 - 53936 p 2 + 6608 p 5 - 32240 p 4 + 22496 p) x 2 
+ (23532 p 3 - 1716 - 12072 p 4 + 2288 p 5 - 21240 p 2 + 9208 p) x 

+ 744 p - 54 + 2812 p 3 - 1552 p 4 + 300 p 5 - 2182 p 2 

I 



and 

v(p) = (1 -p) (2-p) (43 p 3 - 123p 2 + 90p- 27) . 

Although when n = 2 closed form expressions could be 
written down for xo and yo (and hence /oo.2(p)) as a 
function of p, the long expressions are not very infor- 
mative. We can, however, obtain highly accurate nu- 
merical results from Eqs. (fTTj) and (fT2|) for the n = 2 
cylinder lattice for different values of p. For example, 
when p = 1/2, we can solve xo and yo numerically from 
Eqs. |n|) and CE2) as x = 0.389620618156217959 and 
y = 0.442004100446556690, which lead to fooMl) = 
0.643863506776659088. 

For the n = 3 cylinder lattice, in order to solve Xo (or 
yo), we have to solve a polynomial equation with a de- 



gree of 10. When p is not so close to 1, reliable numerical 
solutions can be obtained. For example, when p = 1/2, 
x and y can be solved as 0.441361340073863149 and 
0.277272018269763844 respectively, leading to /oo,3(|) = 
0.632058256526951594. These "exact" numerical values 
can be used to check the results obtained previously by 
the computational methods (Table I, Ref. |2|). This pro- 
cedure confirms the conclusion that the computational 
methods used previously give results with up to 12 and 
13 correct digits. 

For bigger n, it becomes increasingly difficult to solve 
the system of polynomial equations Eq. ([8]). Even nu- 
merical solutions become highly unstable, especially at 
high dimer density. In the following we investigate the 
series expansions of the free energy for lattice strips n x oo 
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for small values of n. Since the behaviors of the solution 
Xq and uq, and hence the asymptotics of the free energy 
foo,n(p), are quite different at the high and the low dimer 
density limits, we discuss the two cases separately. 



III. ASYMPTOTICS AT THE HIGH DIMER 
DENSITY LIMIT 



For clarity we define u = 1 — p. At the high dimer 
density limit when u — > 0, numerical calculations show 
that for both odd and even n, xo — > oo and yo — * 0. If 
we expand l/x and yo as series of u, from Eq. © it is 
found that l/x and y have different leading terms in 
the series expansion for odd and even n. For odd n, the 
leading term of 1 /xq is u 2 and the leading term of yo is 
U n . For even n, the leading terms of X/xq and yo are u 
and u n l 2 respectively: 
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when n is odd 



when n is even. 



These differences in the leading terms of the series expan- 
sions of xo and yo lead directly to the different coefficients 
of u In u in foc,n(u) for odd and even n. By using Eq. JTT 
we obtain for odd n, 
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The difference in the coefficients of win ii in foo,n(u) 
comes directly from the different leading terms of 1/xq 
for odd and even n. 

Some explicit expressions of foo,n(p) f° r cylinder lat- 
tices and lattices with free boundaries are listed in the 
following. 



A. Cylinder lattices 

The asymptotic expansions of foo,n(p) m cylinder lat- 
tices at the high dimer density are listed below for 
n = 1,...,5. For cylinder lattices, the constant term 
of foo,n(p) is given by the following exact expression [|| 
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which can be used to check the constant terms in the 
following results. 
For n = 1, 
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B. Lattices with free boundaries 

The asymptotic expansions of /oo,n (p) in lattices with 
free boundaries at the high dimcr density are listed below 
for n — 1, ... ,4. For lattices with free boundaries, the 
constant term of foo,n(p) is given by the following exact 
expression Q 
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which can be used to check the constant terms of the 
following results. 
For n = 2, 
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ASYMPTOTICS AT THE LOW DIMER 
DENSITY LIMIT 
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Unlike the high dimer density case, at low dimer den- 
sity when p — > 0, numerical calculations show that Xo 
approaches zero and yo approaches 1, for both odd and 
even values of n. The series expansions of Xg and yo thus 
have the following forms: 



From Eq. ([10)) , the general form of the free energy at the 
low dimer density is 



For both odd and even values of n, at low dimer density 
the coefficient of logarithmic term plnp is —1/2, consis- 
tent with the previous results obtained by computational 
methods H . This coefficient comes directly from the fact 
that the leading term in the series expansion of xq is a\p. 

The asymptotic expansions of free energy in cylinder 
lattices at low dimcr density show an interesting prop- 
erty: for lattice strip n x 00 with a width of n, the 
first n terms in the series expansion of /oo,n(p) is exactly 
the same as the first n terms in the series expansion of 
/oo,oo(p)i the free energy of the infinite lattice. In order 
to compare the free energy in semi-infinite 11x00 lattice 
strips with that of an infinite 00 x 00 lattice, in the fol- 
lowing section the series of Gaunt [l(| is used to derive 
the series of /oo,oo(p)- 



A. Free energy for an infinite lattice 

Gaunt gave a series expansion of the dimer activity x 
as a function of the number density t (Ref. IToL column 2 
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of Table II): 

x(t) =t + 7 t 2 + 40 t 3 + 206 t 4 + 1000 < 5 
+ 4678 1 6 + 21336 1 7 + 95514 i 8 
+ 421472 1 9 + 1838680 t w + 7947692 1 11 
+ 34097202 t 12 + 145387044 t 13 
+ 616771148 t 14 + 2605407492 i 15 + ■ • • . 



(13) 



Here t = 6 /4, where 8 is the average number of sites 
covered by dimers when grand canonical ensembles are 
considered @. 

By using the relation between the grand canonical en- 
semble and the canonical ensemble Q and identifying 8 
with the dimer density p, we have the following relation 
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From this relation and the series in Eq. (|13p we can 
obtain the series expression for the free energy of the 
infinite lattice 
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B. Cylinder lattices 



The coefficients of p % in the series expansion of /oo,n(p) 
at low dimer density is listed in Table [J for n — 1, . . . , 7. 
The term — ipln(p), which is common to lattices of all 
sizes, is not included in the Table. Also listed in the last 
row of the Table are the coefficients for the infinite lattice 
(Eq. ([U])). It is evident from the Table that for lattice 
strip nxoo with a width of n, the first n terms of foo,n{p) 
(including the term of — ipln(p)) is exactly the same 
as the first n terms of the infinite lattice /oo,oo(p)- F° r 
example, the lattice strip 7 x oo has the first seven terms 
identical to those of the infinite lattice, up to the term 
of p e . This nice property gives a quantitative estimate 
of the error when we use the values of finite lattices to 
approximate the properties of the infinite lattice. It also 
explains why the sequence of the free energy in cylinder 
lattices converges so fast, especially when p is small 0. 

The term of p n in foo,n(p) is the first term that dif- 
fers from the series expansion of /oo,oo(p). The difference 



between the coefficients of p n in foo,n(p) and foo,oa(p) 
shows a regular pattern: starting from n = 2, the differ- 
ences are 32 , ~ Y92 i 1024 > — 5120 > 24576 ' ~~ 1 14688 ' ^ OT 

example, for n = 2, — 1| — (— j^) — The closed form 
of this alternating sequences clearly is ^— js-- 

The difference between the coefficients of p n+1 , the sec- 
ond term that differs between finite and infinite lattices, 
also shows a regular pattern. Starting from n — 3, the 

sequence is — , JQ2Z > — 4096" > 16384 » — 65536 > " ^ s 

vious that the sequence takes the closed form expression 

^Vv n ■ Due to the limited number of data points in Ta- 
ble HI currently it is not clear whether the differences of 
higher degree terms also can be written in simple closed 
forms. 

From the closed form expressions of these two domi- 
nant terms that differ between finite and infinite lattices 
at low dimer density, we see that when n > 3 



/oo,oo(p) ^ foo,n(p) 
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+ 0(p n+z ). 

The coefficients of p 1 in the series expansion of xq and 
— (lnyo)/n are also listed in Tables HT1 and fTTTl From the 
Tables we can see that for cylinder lattice strips nxoo, xq 
and (hiyo)/n share the same first n — 1 terms with their 
corresponding series expansions of the infinite lattice. 



C. Lattices with free boundaries 

The coefficients of p l in the series expansion of foa,n{p) 
for lattices with free boundaries at low dimer density 
is listed in Table IIVI for n = 2,3, and 4. The term 
— ^pln(p), which is common to lattices of all sizes, is 
not included in the Table. From the Table we can see 
that, unlike cylinder lattices, none of the coefficients is 
the same as the coefficients of the infinite lattice. 

As for the cylinder lattices, the coefficients of p l in the 
series expansion of xo and —(\nyo)/n are also listed in 
Tables |V] and IVI1 Again, except for the coefficient of p 
in (lnj/o)/ n j none of the coefficients is the same as the 
coefficients of the infinite lattice. 



APPENDIX A: BIVARIATE GENERATING 
FUNCTIONS OF MONOMER-DIMER MODELS 
IN TWO-DIMENSIONAL PLANAR LATTICES 

The bivariate generating functions of monomer-dimer 
models can be derived directly from the matrices M n 
which are used in the computational studies @, 0, H, [1] • 
The denominator H(x, y) of the bivariate generating 
functions G(x, y) is directly related to the characteristic 
function of M n : 

H(x,y)=det(M n -y- 1 I)xy w , 

where w is the size of the matrix M n . For the purpose 
of this article, it is sufficient to know the information 
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TABLE I: Coefficients of p' in the series expansion of /oo,n(p) of cylinder lattices at low dimer density. The term — |pln(p), 
common to lattices of all sizes, is not shown. The last row for the infinite lattice /oo,oo(p) is taken from the series expansion 
Eq. (|14|l . Terms of /oo,n(p) that are equal to those of the infinite lattice are underlined. 



n 


P 


-P 2 


-P 3 


-P 4 


-P 5 


-P 6 


-P 7 


-P S 


1 


ln2+l 
2 


3/8 


7/48 


5/64 


31/640 


21/640 


127/5376 


255/14336 


2 


In 2 + 1/2 


13/32 


115/768 


419/6144 


5491/163840 


17489/983040 


116687/11010048 


423771/58720256 


3 


In 2 + 1/2 


7/16 


1/6 


127/1536 


1027/20480 


8653/245760 


1 7(177 /fifift 1 OS 


9^81 /l '"tl (179 


4 


In 2 + 1/2 


7/16 


31/192 


239/3072 


461/10240 


3569/122880 


27325/1376256 


7579/524288 


5 


In 2 + 1/2 


7/16 


31/192 


121/1536 


473/10240 


941/30720 


3791/172032 


4375/262144 


6 


In 2 + 1/2 


7/16 


31/192 


121/1536 


471/10240 


1243/40960 


3707/172032 


4177/262144 


7 


ln2 + 1/2 


7/16 


31/192 


121/1536 


471/10240 


1867/61440 


3719/172032 


4215/262144 


oo 


In 2 + 1/2 


7/16 


31/192 


121/1536 


471/10240 


1867/61440 


7435/344064 


4211/262144 


TABLE II: The coefficients 


in the series expansion of xo 


for cylinder lattice strips n x 


oo at low dimer density. 


Terms that are 


equal to those of the infinite lattice are 


underlined. 










n 


P 


P 2 


P 3 


P 4 


P 5 


P 6 


P 7 


P 8 


1 


1/2 


3/4 


l 


5/4 


3/2 


7/4 


2 


9/4 


2 


1/4 


13/32 


71/128 


1393/2048 


6353/8192 


55073/65536 


230343/262144 


7519577/8388608 


3 


TJT 


7/16 


81/128 


423/512 


4179/4096 


9993/8192 


23341/16384 


854147/524288 


4 


1/4 


7/16 


5/8 


411/512 


497/512 


289/256 


20917/16384 


370861/262144 


5 


1/4 


7/16 


5/8 


103/128 


2001/2048 


9369/8192 


42803/32768 


192145/131072 


6 


1/4 


7/16 


5/8 


103/128 


125/128 


9355/8192 


21329/16384 


190871/131072 


7 


1/4 


7/16 


5/8 


103/128 


125/128 


2339/2048 


42673/32768 


191043/131072 


oo 


1/4 


7/16 


5/8 


103/128 


125/128 


2339/2048 






TABLE III: The coefficients 


in the series expansion of — 


ln(y )/n for 


cylinder lattice strips n x oo at low dimer density. Terms 


that 


are equal to those of the infinite lattice are underlined. 








n 


p 


P 2 


P 3 


P 4 


p 5 


P 6 


P 7 


P 8 


1 


1/2 


3/8 


7/24 


15/64 


31/160 


21/128 


127/896 


255/2048 


2 


TJT 


13/32 


115/384 


419/2048 


5491/40960 


17489/196608 


116687/1835008 


423771/8388608 


3 


TJT 


7/16 


1/3 


127/512 


1027/5120 


8653/49152 


17677/114688 


16667/131072 


4 


1/2 


7/16 


31/96 


239/1024 


461/2560 


3569/24576 


27325/229376 


53053/524288 


5 


1/2 


7/16 


31/96 


121/512 


473/2560 


941/6144 


3791/28672 


30625/262144 


6 


1/2 


7/16 


31/96 


121/512 


471/2560 


1243/8192 


3707/28672 


29239/262144 


7 


1/2 


7/16 


31/96 


121/512 


471/2560 


1867/12288 


3719/28672 


29505/262144 


oo 


1/2 


7/16 


31/96 


121/512 


471/2560 


1867/12288 







TABLE IV: Coefficients of p 1 in the series expansion of /do, n(p) of lattices with free boundaries at low dimer density. The term 
— ^pln(p), common to lattices of all sizes, is not shown. Numbers in square brackets denote powers of 10. 







p 


-P 2 


-p 3 


-p 4 


-P 5 


-p 6 


-p 7 


-P 8 


2 




(Jn 3+ 1) 


5/12 


17/108 


49/648 


403/9720 


25/972 


3223/183708 


27569/2204496 


3 




,(10/3) + 1) 


S7/2[2] 


1569/1[4] 


14697/2[5] 


2208627/5[7] 


192832569/625[7] 


9382406661/4375(8] 


1278517726503/875[11] 


4 


Jo 


n(7/2) + 1) 


43/98 


1123/7203 


26323/352947 


187546/4117715 


18361592/605304105 


850524016/41523861603 


70713460952/4747561509943 



TABLE V: xo The coefficients in the series expansion of xo for n x oo lattices with free boundaries at low dimer density. 
Numbers in square brackets denote powers of 10. 



n p 


P 2 


p 3 


P 4 


P° 


p 6 


p 7 


P 8 


2 1/3 


5/9 


7/9 


239/243 


851/729 


2909/2187 


29017/19683 


31507/19683 


3 3/10 


261/5[2] 


9207/125[2] 


116397/125[3] 


4826031/3125[4] 


021574591/15625[5] 


28645564383/1953125[4] 


3178995484149/1953125[6] 


4 2/7 


172/343 


11888/16807 


738720/823543 


6219872/5764801 


354308800/282475249 


137557926784/96889010407 


7481262371584/4747561509943 



TABLE VI: — ln(yo)/n The coefficients in the series expansion of — (\nyo)/n for n x oo lattices with free boundaries at low 
dimer density. Numbers in square brackets denote powers of 10. 





P 


P 2 


P 3 


P 4 


p 5 


p 6 


P 7 


p 8 


2 


1/2 


5/12 


17/54 


49/216 


403/2430 


125/972 


3223/30618 


27569/314928 


3 


1/2 


87/2[2] 


1569/5(3] 


44091/2[5] 


2208627/125[5] 


192832569/125[7] 


28147219983/21875 [7] 


12785 17726503/ 125 [11] 


4 


1/2 


43/98 


2246/7203 


26323/117649 


750184/4117715 


18361592/121060821 


1701048032/13841287201 


70713460952/678223072849 
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of H(x,y). From H(x,y) and a few initial terms, the 
bivariate generating functions can be obtained. For com- 
pleteness we list G(x, y) for small values of n in this Ap- 
pendix. The cylinder lattices and lattices with free edges 
are listed separately. Variable x is associated with the 
number of dimers s and variable y is associated with the 
length of the lattice m, as defined in Eq. ([7]). 

When n = l, there is no distinction between these two 
boundary conditions: 



1 — y — xy 2 



1. Cylinder lattices 

For n = 2, 

G c = i-yx 

2 1 — (1 + 3 x) y + x (x — 1) y 2 + y 3 x 3 

For n = 3, 



G c = l-2yx- y 2 x 3 

3 x 6 y 4 - x 3 (x - 1) y 3 - x (1 + 3 x + 5 x 2 ) y 2 - (5 x + 1) y + 1 ' 



= x 6 y 4 + x 4 y 3 - 2 x 2 (1 + x) y 2 - xy + 1, 



I 

For n = 4, G\ = F 4 /H 4 , where For n = 3, = F^/Hf, where 

F 4 = -x 8 y 4 + 3x 5 y 3 + 4x 4 y 2 - x(3 + 4x)y + l, 

and 

H A =x 12 y 6 - x 8 (-x + 2x 2 + l)y 5 and 
-x 5 (6x 2 + 2x + l + 9x 3 )y 4 

+ 2x 3 (l3x 2 + 5x + 1 + Ax 3 ) y 3 flf = - x 9 y 6 + x 6 (x - 1) y 5 

+ x (-6 x - 1 + 7 x 3 - 6 x 2 ) y 2 + x 4 (5 x 2 + 3 x + 2) y 4 + x 2 (2 x + 1) (x - 1) y 3 

- (x + 1) (6 x + 1) y + 1. - a; (1 + x) (5 as + 2) y 2 - (1 + 3 x) y + 1. 

For n — 5, G§ = F 5 /H 5 , where 

F 5 = - x 15 / + 2," (-2 + x) y* For n = 4 ' ° 4 = with 

+ x 8 (8x 2 + 2 + llx)y 4 rfb i4„,7 ii..e ^8/. ,,2 



+ 2x 5 (7x 2 + 3 + 8x) y 3 
-x 3 (2-x + 8x 2 )y 2 
-2x(2 + 5x)y + l, 



and 



H 5 =x 20 y s + x 15 (3x 2 -x + l)y 7 



iPp =x 14 y 7 - 2 x 1 V - a; 8 (3 x + 5 x 2 + 3) y 5 
+ x 6 (2x + 3)(l + x)y 4 
+ x 4 (2x + 3)(l + 3x)y 3 
-3x 2 (l + 3x + x 2 )y 2 
-2x(l + x)y + l, 



x 11 (l9x 4 + llx 3 + 7x 2 + 2x + l)y 6 and 
x 8 (2 x 4 + 65 x 3 + 39 x 2 + 1 1 x + 2) y 5 

x 5 (41 + 95 x 4 + 39 - 9 x 2 - 6 x - l) y 4 < = - x 18 y 9 + ^ W H + ^ + 4 ) V 8 

x 3 



(34x 4 + 85x 3 + 69x 2 + 19x + 2) y 3 + x 11 (2 + 9 x 3 + 3 x + 9 x 2 ) y 7 

■4 , in,3, 27 J, 1n - +l N, J 3 -x 8 (l9x 3 -l + 6x 2 + 5x 4 )/ 



a; (19 x 4 + 19 x 3 + 27 x 2 + 10 as + l) y 2 

(l5x 2 + 9x + l)y+l. -x 6 (29x 2 + 3 + 14x + 24x 3 + 21x 4 )y 5 

+ x 4 (41 x 2 + 40 a; 3 + 9.x 4 + 18 a; + 3) y 4 



Lattices with free boundaries 



x 2 (4x 2 - 1 - 4x + 15a; 4 + 27 x 3 ) y 3 



- x (5 a; 3 + 2 + 21 a; 2 + 13 x) y 2 
For 7i = 2, + (-1 - 3x 2 - 5x) y + 1. 



G fb _ l - x y 



x 3 y 3 — xy 2 — (2 x + 1) y + 1 
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